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Description 

BACKGROUND OF THE INVENTION 

The present invention relates to a method of producing a train running plan which is effective for realizing a high 
speed, high density, and energy-saving running of a railway train. 

With respect to a method of producing the optimal train running plan under the conditions such as line conditions, 
e.g., a limit speed and a grade, and rolling stock characteristics, e.g., a running resistance, driving fbrce characteristics, 
and speed decreasing cheiracteristics. there are puk)lished the following literatures: 

Literature 1 ; 5th National Symposium for Utilization of Cybernetics in a Railway, June, 1968. pp. 11-16: and 
Uterature 2; The Transactions of the Institute of Electrical Engineering Engineer's o1 Japan. VoL 106-B. No. 9, Sep- 
tember. 1986, ppi 769-776. 

The method in the literature 1 is such that the optimal running method is obtained in the running of arbitrary form 
between stations on the basis of the dynamic planning by treating the running time, the power consumption and the 
number of notch switching as the objective fectors. This method is obtained in the form of notch sequence. Then, the 
running time, the power consumption and the like as the objective factors are obtained by performing the running sim- 
ulation of tiie train to determine the position, tiie speed and the like of the train. 

In the literature 2. the discussion is made with respect to the energy-saving running method in the running of the 
super-express rallway*s train stopping every station, between the stations. In the literature 1 . the calculation is made by 
the simulation with the time base being divided since the equation of motion of the train cannot be exactly obtained in 
the leu'ge area. On the other hand, in the literature 2. each of the strict solutions is obtained on the assumption that the 
magnitude of the grade takes a fixed value and the resultant solutions are linked to one another to perform the calcula- 
tion of ttie consumed energy, the necessary time and the like. 

The discussion about the energy-saving running method in the literature 2 is made in such a way that four patterns 
near the actual running are prepared, the consumed energy is calculated every running pattern on the assumption that 
the magnitude of the grade Is zero over the whole territories, and the predetermined comparison is performed. Out of 
the four running patterns, the running pattern having the minimum consumed energy provkies the results of the maxi- 
mum acceleration -> the fixed running at the maximum speed -> stopping by the repetition of the coasting and the nor- 
mal maximum kxaking. 

in the alxxve-mentioned running pattern -having the minimum oonsjmed energy, only the maximum speed is 
changed to obtain the speed in whk^h the necessary time coir^es with the schedule time for realizing the scheduled 
running in which the necessary time coincides with the schedule time. 

(However, in the prior art method of obtaining the optimal running plan on the basis of the dynamic planning in the 
literature 1 , there arises a problem in that the calculation requires much time. 

In the literature 2. as the running pattern of the Gmit speed between ttie stations, only the pattern shown in Rg. 1 is 
supposed. Thus, the method described tiierein can be applied to only this running pattern. Therefore, it cannot be 
applied to the complicated pattern as shown in Fig. 2 for example. 

Further, in addition to the literatures 1 and 2, with respect to the method of changing the train running plan pro- 
duced in advance in the train in a real-time manner, there is published the literature 3. i.e.. The Transactions of ttie Insti- 
tute of Electrical Engineering Engineer's of Japan. Vol. 107-D, No 5. May, 1987. pp. 665-672. 

The method described in the literature 3 is about the runrting plan of the super-express railwa/s train stopping 
every station. In this case, the running pattern as shiown in Rg. 1 is supposed as the running pattern of the limit speed. 
Then, out of a plurality of running methods which are supposed in advance, one providing the minimum energy-saving 
is determined on the assumption that the magnitude of the grade is zero over the territory shown in Fig. 1 . The resultant 
energy-saving running method is such that the running pattern is obtained by combining the maxinuim acceleration, the 
fixed speed running, the coasting and the maximum deceleration with one another. 

In the running method having such a running pattern, the parameters which are to be determined in advance are 
the maxinnum speed and the deceleration starting point in the running territory. The two parameters are changed so as 
to provide the maximum energy-saving running under tine condition of the scheduled running (the running time coin- 
ckfes with the schedule time). 

In tiie method described in the literature 3. only the running pattern as shown in Rg. 1 is supposed as the running 
pattern of the limit speed between the stations, the method of energy-saving running isdeternrtined in advance, arxi only 
the determination of the deceleration starting point is performed in the real-time chiange in tfie train. Therefore, this 
method cannot cope with the occurrence of change of the Dmit speed as shown in Rg. 3A. for example, which is not 
supposed in advance. 
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SUMMARY OF THE INVENTION 

It is the object of the present invention to determine the target speed so as to provide the minimum consunned 
energy during the running of the train. 

This otaject is met by a method according to claim 1 . Preferred embodiments are disclosed in the depending claims. 

The general pattern of the limit speed between the stations is supposed as shown in Rg. 4, and the running pattern 
on the basis of that limit ^eed is assumed here. Then, the limit speed territories in Fig. 4 are numbered serially (the 
number of tenrltorles is N). At the same time. It is assumed that the limit speed in the l-th limit speed tem'tory is V^^x i 
and the running speed to be targeted (hereinafter, refen'ed to as simply Ihe target speed", when applicable) is VI. Then, 
the relationship of Vi ^ ^max,\ established. 

Under ^ch assumptions, to set the target speed (it is not always possible to reach this target speed) every limit 
speed territory is determined the train running method between the associated stations. For example, in the running 
pattern of the maximum acceleration -> the fixed speed running -> the deceleration by the normal maximum braking, 
under the limit speed of Vmaxj 0=1 N). if Vi (1=1 N) is established, the train running method between the asso- 
ciated stations is correspondingly established. Conversely, it is also true that to establish the train running method 
between the stations is to establish Vi (i=1. N). 

By the optimal running method in the train running method as described above, it means the scheduled running in 
which the necessary time required for the train to run between the stations coincides with the schedule time and at this 
time, the minimum consumed energy is obtained. 

It Is assumed that the necessary time required for the train to run through a predetermined territory Is T and the 
consumed energy Is E. Then, the necessary time T and the consumed energy E cannot be analytically calculated even 
if the line conditions and the rolling stock characteristics as already described are known. However, these factors can 
be obtained by making the equation of motion of the train discrete corresponding to the travel distance and by perform- 
ing the numerical calculation (i.e., by performing the simulation). 

At this time, if Vi (1=1 N) is established, the necessary time T and the consumed energy E are uniquely obtained. 

Therefore, each of T and E is a function of Vi (i=1 N). That is. the following relational expressions are established: 

T = T(Vi Vn).E = E(V, V^) (1) 

Then, to obtain the optimal Vr (i=1 , ... , N) is established mathematically in such a way as to obtain a set of target speeds 

Vi, Vn for minimizing the objective function E(Vv .... V^). under the condition of the necessary time T(Vi V|sj) = 

the schedule time Tq. 

More specifically, Vi is changed so as to decrease the consumed energy E, and for example, the partial, derivative 
8E/dV\ of E witii respect to Vi is obtained. Thus, the direction of decreasing of E is obtained. There is known a non-linear 
planning in which the calculation for obtaining the direction of decreasing of E is repeated until the necessary time T 
approaches tiie scheduled time Tq within ttie range of a certain allowable error 6T, and when the necessary time T 
becomes in the range of the allowable error, tiie above calculation is finished. 

Or. the following equation may be used as the ot^ective function. 

J«E + a{(T-TD)/5T} (2) 

In this case, it is possible to simultaneously evaluate both the consumed energy E and the necessary time T. At this 
time, the evaluation by only the consumed energy E can be performed in the case of a » o. 

The above description has been given with respect to the specific case where tiie train runs through the territory 
between tiie two stations, arvd it is then stopped at the following station. However, even in tiie case where both the two 
stations are passed stations, the same means can be used by giving the initial speed at the first station in addition to 
the above conditions. 

When the mnning speed of the train is changed, tiie subsequent target speed y » (V^, V2 V^) as a vector form 

is newly obtained on tiie basis of tiie change ratio of tiie consumed energy between that running speed and tiie target 
speed of interest. Therefore, tiie component of the gradient vector vE(y) in the direction of decreasing of the consumed 
energy E is obtained from the following approximate expression. 

aE/aV(n) ^ AE/AV(n) (3) 

That is, the target speed V(n) is changed in the direction of tiie gradient vector, and tiie target speed V(n) for minimizing 
the consumed energy E is obtained. 

According to another embodiment of tiie present Invention, Instead of ttie gradient vector vE(y}, tiie ratio of the 
decreasing of the consumed energy E to the increasing of the running time T when V(n) is changed by AV(n). i.e.. the 
ratio of AE to AT is used. Then, each of the components of each of tiie gradient vectors is obtained from tfie following 



3 



EP0 467377B1 

equation. 

-AE/AT = -{AE/AV(n)}/{AT/AV(n)} (4) 

= -{dE/avcnMaT/avtn)} 

5 

The value of the equation (4) is called p(m). Then, let p(n) be replaced with a as the form of vector. Then, the target 
speed y Is changed in the direction of q, In this case, it is assumed that the initial value of the target speed V(n) is 
expressed by {the limit speed Vmax(n) - (the marginal speed)}. Then, by the marginal speed it means the upper limit 
speed that does not exceed the limit speed during the running at the fixed speed. 

10 In the calculation of the running time, in addition to the running curve on the speed curve which is expressed by the 
speed on the axis of ordinate with respect to the position on the axis of abscissa as shovm In Rg. 5. the curve of the 
inverse number of the speed is also utilized. Since the single integral of tiie Inverse number of tiie speed witii respect 
to the position on the axis of abscissa becomes the running time of the territory of interest, it is herein referred to as the 
running time density", and is expressed by A. 

IS The characteristics of the train are characterized by the acceleration and deceleration characteristics, the running 
resistance and tfie like. With the train, it is assumed that all the grades of the traffic line are zero. Then, as shown in Rg. 
7. in the acceleration up to the suitable speed more than the maximum speed on the line, the data about the running 
speed y. the running time density X and the running time T with respect to the running disteince x are produced. Now. 
the running speed v. the running time density X and the running time T are expressed by the following equations, 

20 respectively. 

v = p(x) (5) 
>. = 1/p(x) = f(x) (6) 



2ff 



30 



X 



40 



T-Jf(Odi;«F(x) (7) 

0 

Moreover, as shown in Rg. 8. the individual data in the deceleration from the running speed v of interest are also pro- 
duced. In this case, the running speed v, the running time density K and the running time T are expressed by the follow- 
ing equations, respectively. 

v = q(x) (8) 
X = Mqix) = g(x) (9) 

X 

T=Jg(C)dC = G(x) (10) 



Since each of the functions p(x). q(x), f(x). g(x). F{x) and G(x) is a monotonous function, these functions have 
45 respective inverse functions. Moreover, as the furtctions necessary for the calculation of tiie running time, the data of 
F(f ^ (X)) with respect to the variable X in the acceferation and tiie data of G(g'^ (X)) with respect to the variable X in the 
deceleration are produced. 

Although these functions are not always dearly expressed in the form of formula, they are stored in the computer 
in tiie form of discrete numeral data. 

so Fig. 9 is a graphical representation showing the runrting of the train in which tiie fixed speed running is performed 
up to the position x^ at tiie speed . tiie acceleration is started from tiie position x^. tiie running speed reaches \he 
speed V2 in the positfon . the fixed speed running is performed up to the position X2 at the speed V2. the deceleration 
is started from the position X2, the running speed reaches the speed V3 in the position ^2. and the fixed speed running 
is perfomied up to tiie position X3. Incidentally, in ttiis case, it is also assumed that the magrvtude of tiie grade is zero 

55 in each of the running territories. 

Moreover, Rg. 1 0 is a graphical representation showing the inverse nunrfoer of the speed. i.a, the running time den- 
sity. The running time density is expressed by the following equation. 

Xi=1/Vi(i=1.2.3) (11) 
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The position where the running is switched over from the acceleration to the fixed speed running can be obtained 
using the speed curve p(x) in the acceleration of Rg. 7, and the distance required for the acceleration from Vi to Mz can 
be obtained by the following equation. 

s P'V2)-pV,) = ^,-x, (12) 

\- The running time T^ from to %^ corresponds to the area A shown in Fig. 10. In order to obtain this running time 
Ta; the travel territory tsetween xi and %2 converted into the characteristics shown In Fig. 7 and the resultant charac- 
teristic curve of Fig. 7 is then integrated. Rrst the positions y^ and y2 corresponding to X\ and X2 are given by 
10 y 1 s f ^ } and y 2 f \^ 2) > respectively. Therefore, the running time Ta is given by: 



" ' TA = Jf(C)dC= J f(OdC 

IS 

= F{f*\>.2))-Fa''{>-i)) 



Moreover, the running time Tb from to X2 corresponds to the area B and is given byX2*(x2-^i).Asa result, 
20 the running time Tp^ from x^ to X2 is obtained using the following equation. 

Tab = Fff '0.2)) - FP 'O^ 1)) + ^2 • (^2 ' ^1) (14) 

The position ^ where the running is switched over from the deceleration to the fixed speed running is obtained firom 
25 the followving equation, using the speed curve q(x) in the deceleration of Fig. 8. 

q ^V3)-q-V2) = 42-X2 (15) 

in the same manner as in the acceleration, the running time Tc from to Is given by: 

30 

f g(OdC 

^ =G(g\x3))-G(g\x2)) 

The running time Tq from %2 ^ X3 conresponds to the area D and is given by * (X3 - 4 2) • ^ ^ result the run- 
ning time Tqo from X2 to X3 is obtained using the following equation. 

40 

T CD = G(9 (>. 3)) - G(g ■\>^2)) + ^ 3 • 3 - ^ 2) (17) 

The resison for performing tiie calculation on the assumption ttiat the magnitude of the grade is zero in the above 
description is that since the magnitude of the grade depends on the place, taking the actual grades into consideratton, 
45 it is impossible to produce the desired data in advance by the simulation. However, tiie following items are true: 

1 . Since the recalculation of the running time in the train is necessary for the territory where the change off the limit 
speed occurs locally and tiie subsequent several territories, the distance requiring the recalculation is short. Thus, 
there is a snrrall error between this case and the case where the magnitude of the grade is taken into consideration. 
so 2. Such a territory where the change of tiie limit speed occurs is mainly included in the high speed area. arxJ tiiere- 
fore, tiiere is sinalt influence of some variation of the speed due to the grade upon the running time. 

For tiie above reasons, there is tiie small error between the present case and the case where the magnitude of the 
grade is taken into consideration. 
55 As described above, the necessary data are calculated in advance by the simulation and the like, and the running 
time is then calculated using the resultant data, whereby it is possible to correct the running method which is set in 
advance, in a real-time manner. 

According to the present invention, the necessary data is calculated in advance by the simulation Even when the 
necessity of correcting the train running method occurs, e.g., the change of the limit speed temporarily occurs, it is pes- 
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sible to oorrect the running method in the train in a real-time manner. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 Rg. 1 is a graphical representation showing a prior art running pattern of a limit speed between stations; 

Fig. 2 is a graphical representation showing an example of a running pattern of a limit speed between stations; 
-Rg. 3A is a graphical representation showing an example of the running pattern in the case where the setting of 
the limit speed occurs temporarily; 

Rg. SB is a graphical representation showing an example of the running pattern con^esponding to the inverse 
10 number of the limit speed illustrated in Rg. 3A; 

Rg. 4 is a graphical representation showing an example of the running pattern of the limit speed between the sta* 
tions; 

Rg. 5 is a graphical representation showing a running pattern used for correcting the running plan; 
Rg. 6 is a graphical representation showing the inverse number of the running pattern shown in Rg. 5; 
IS Fig. 7 is a graphical representation showing the speed, the inverse number of the speed and the running time in the 
acceleration in the case where the magnitude of the grade of the travel tenitory is "0"; 

Rg. 8 is a graphical representation showing the speed, the inverse number of the speed and the running time in the 
dec^eration in the case where the magnitude of the grade of the travel territory is "0"; 

Rg. 9 and Fig. 10 are graphical representations useful in explaining the calculation of the running time using the 
20 Inverse number of the speed; 

Fig. 1 1 is a biockdiagram showing the an^angement of a system forcarrying out a method of producing a train run- 
ning plan according to the present invention; 

Rg. 12 is a flow chart showirtg the method of producing a train running plan of a first embodiment according to the 
present invention; 

2S Rg: 1 3 is a diagram useful in explaining an example of the running plan which is displayed on a CRT of a train sys- 
tem; 

=Rg. 14 is a graphical representation showing an example of the grade data of a limit speed territory; 
Fig. 15 is a graphical representation showing the running pattern of the limit speed which is used when a train 
passes a station; 

30 Rg. 1 6 Is a flow chart showing the method of producing a train running plan of a second embedment according to 
the present invention ; 

Rg. 1 7 is a flow chart showing the method off correcting a train running plan; 

Fig. 18 and Rg. 19 are graphical representations useful in explaining examples of correction of the running plan; 
and 

35 Rg. 20 and Rg. 21 are graphical representations useful in explaining another examples of correction of the running 
plan. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

40 Preferred embodiments of the present invention will hereinafter be described in detail with reference to the accom- 
panying drawings. 

Rg. 1 1 is a biockdiagram showing the arrangement of a system for carrying out a method of producing a train run- 
ning plan according to the present invention. 

This system arrangement is made up of a train system 200a having a main function of displaying an optimal train 
45 running plan on a C RT 203 provided in a train to perform the operational support for an engineer, and a ground system 
200b having a main function of producing the optimal train running plan to be displayed by tiie train system 200a. 

Moreover, for brevity, it is assumed that the traffic route is a line consisting of two rails, a plurality of stations are 
included in the ti^ffic route, a running train is limited in type, and each of the trains is a local t^n. 

A running plan producing unit 208 serves to refer to ttie necessary data in a rolling stock characteristic data file 21 2 
so and a line condition data file 21 3 to produce optimal train running plan thereby to store the optimal train running plan in 
a running plan data file 21 1 . This processing is performed in advance witii respect to all the territories, the type of the 
train, the up train line and the down train line. 

Then, to determine the train plan is to determine a target speed in each of the limit speed territori^ Incidentally, in 
addition to the target speed, the curve of the speed corresponding to tiie position on the line, which is designated ttie 
55 running curve, will be also displayed in tiie train. In the data file 21 1 , tiie result of the simulation when the optimal target 
speed is set will be recorded in the form of the position and the speed. The details of the method of producing an optimal 
train running plan will be described later with reference to Rg. 12. 

Then, by the rolling stock characteristic data it means the drive force characteristic data and the deceleration char- 
acteristic data corresponding to the train speed, the running resistance which is a function of the second order of tiie 
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train speed, the train organization, the train weight and the like all of which are inherent in the train. These rolling stock 
characteristic data are stored corresponding to the type of the running train. 

Moreover, by the line condition data it means the position of the station (stopping target position), the limit speed 
information (territory starting position and limit speed, eta), the grade information (starting position and magnitude of 
5 grade, etc.) and the like all of which are obtained with the starting point of the traffic route being treated as a standard. 
Since it is assumed that the traffic route is a single line, these data are limited to only two kinds with respect to the up 
train tine and the down train line. 

Out of the data of the optimal train running plan thus produced, the data oonresponding to the running train and the 
teaitory t}etween the stations at which the train is to be stopped in order need to be transferred to the train system 2O0a 
10 of the running train. As for the means for transferring such data, there are known various kinds of ones. However, par- 
ticularly in the present example, an IC card is used to transfer tiie data from the ground system 200b to the train system 
200a. The running plan producing unit 208 serves to search for tiie necessary data to take out them from the data file 
21 1 . thereby to send them to an IC card writer 21 4. The IC card writer 214 serves to write the data thus sent thereto in 
an IC card 215. The IC card 215 having the data written therein is the same as an IC card 206 in the t^ain system 20Oa. 
IS A running support unit 201 serves to read out the running plan from the IC card 206 through an IC card reader 207 
prior to the start of the train and to record the running plan in a running plan data file 204 in a train computer. Further, 
the running support unit 201 performs the display of the present position, the target speed at the present position, the 
running curve in the running plan, and the like on a CRT 203. corresponding to the travel position of the ti'ain. A key- 
board 202 in the train system 20ba and the 209 in the ground system 200b are used for the input of the commands and 
20 the like when the respective systems are operated and so forth. An example of the displayed picture is shown in Fig. 1 3. 

In Fig. 13. the reference numeral 701 designates the limit speed, the reference numeral 702 designates the optimal 
target speed, and the reference numeral 703 designates ttie running curve which is obtained by performing the simuta- 
tion of the running at the optimal target speed. The reference numeral 704 indicates tiie present position on the running 
curve. The target speed, the present speed, the present position arKi the running time are numerically displayed in 
25 areas 705. 706. 707 and 708, respectively. 

Returning to Fig. 11. the position of the train is calculated In such a way tiiat a rolling stock position detection unit 
205 nieasures the rotational frequency of tiie wheel and the connection is performed in tiie place where tiie positions of 
the station and the vicinity thereof can be accurately grasped. 

Rg. 12 is a flow chart in accordance with which the running plan producing unit 208 in the ground system 2O0b 
30 shown in Fig. 1 1 carries out a first embodiment of a method of producing an optimal train running plan. 

The first embodiment is characterized in that the target speed is determined so as to provide the minimum con- 
sumed energy when the train is running. 

In Step 101 . the data inputted from the keyboard 209 is read out. The data which are inputted from the keyboard 
209 are made up of the type of the train, the identifiers with respect to the individual stations, the starting time, and the 
3S estimated time of arrival. 

In Step 1 02. on tiie basis of the data thus inputted, the line condition data with respect to the stations at which the 
running train is to be stopped in order (flie line iengtii. the station position, tiie limit speed information, ttie grade infor- 
mation and tfie like) and the rolling stock characteristic data of the running train (the drive force characteristics, the 
deceleration characteristics, the train weight and the like) which are stored in advance in the files 213 and 212 of the 
40 ground system 200b. respectively, are read out. The train weight is calculated on the basis of the average seat-load fac- 
tor (it is estimated from the past data and the like). 

In Step 103, on the basis of the above data thus read out, the target speed which is ol>tained by subtracting the 
marginal speed from the limit speed is assigned to each of the limit speed territories. The marginal speed is set so as 
not to exceed the limit speed during ttie fixed speed running. Thus, for example, if the marginal speed is set to 3 km/h. 
45 when the limit speed is 100 km/h, 97 km/h is set to the target speed. 

In Step 104. the necessary time T and the consumed energy E in tiie case where the train will run between the sta- 
tions at the target speed thus set are calculated. The concrete metiiod of calculating tiie necessary time T and the con- 
sumed energy E will hereinbelow be described. 

The running pattern of tiie limit speed shown in Fig. 4 is used as that in the present embodiment. The axis of 
so abscissa designates the distance from tiie starting station and will hereinbelow be represented by x. The limit speed 
territories are numbered serially (the total number of territories is N). At this time, it is assumed that the starting position 
of the i-th territory is Xj and the limit speed thereof is V|viAx,i- "^^t is. if the position of the treun satisfies the relationship 
of X| £ X ^ Xj^.1. the limit speed at that position is Vmaxj- 

Moreover, as described above, It is assumed that the running pattern takes the combination of the meudmum accel- 
55 eration the fixed speed running the deceleration by the normal maximum braking. 

Next, it is assumed that the grade of the n-th limit speed territory is as shown in Rg. 1 4 and the number of territories 
each having the grade is Mn. Tlien. In tiie position satisfying the relationship of ^.^ ^ x s ^ rn> grade is yp Ind- 
dentally. the unit of the grade is expressed by o/oo (per mil). 

The running resistance Rp is a function of the second order of tiie train speed v and is expressed by the following 
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RR(v) = a + bv + cv^ (18) 



where a, b and c are parameters inherent in the kind of the train. 

The grade resistance Rq is a function of the position x. Since the unit of the grade y is given by o/oo, if the angle 
con^esponding to.the grade is given by 6, the following relationship is established. 

sine = e = tane = y/1000 

Thus, if the train weight is given by M, the following relationship is established. 

RG(x)sM*g«sine »M-g«Y/1000 

15 Then, if the drive force characteristics and the deceleration characteristics corresporKling to the speed v are 
expressed by Tq(v) and Tb(v}. respectively, the equation of motion of ttie train motion are given by: 



M-^ = FM{v)-Rn{v).R3(x) (19) 



where FmCv) is tiie ouiput of a motor_(TQ(v) in the acceleration and -Tb(v) in the deceleration). Since the limit speed as 
the constraint condition depends on the position x. the independent variable of tiietime t is inconvenient to treat. Then, 
the^time lis subjected to tiie changaof variables to obtain the independent variable x. At this time, the equation is rewrit- 
es ten as follows: 

"^•v^^FmM-RrM-RqCx) (20) 

30 

Thereafter, the necessary time T arvl the consumed energy E can t>e obtained in such a way that the distance x is 
made to be discrete and the equation (20) is solved by the numerical calculation. 

Incidentally, alttiough it is assumed that with respect to the train running, the acceleration is performed in accord- 
ance with Tq(v) and the deceleration is performed in accordance with Tb(v), the control in the fixed speed running is 
35 performed in the following manner. 

The lower limit AVd and the upper limit ANAj of the error between the present speed y and the target speed Vn are 
set Then, the control in the fixed speed running is performed in accordance with the following relationship. 



V ^ Vn - AVd : acceleration due to Tq(v) 

Vn - AVd < V < Vn + AVu : coasting (F vi(v) « 0) 

Vn AVu s V : deceleration due to T3(v) 

45 In order to otnain the partial derivative dE/aVi of the consumed energy E with respect to the target speed Vi, the 
target speed Vi is reduced by the small speed eV to perform tiie same calculation. If the consumed energy at this time 
is given by E\ the following relationship Is established. 



dE/aVi = (E - EVeV (21) 



Such a calculation is repeatedly performed with respect to i (i=1 . ... N). 

After tiie completion of the atxsve processing, with the target speed being equal to the initial setting value (i.e.. in 
the f irst calculation), tiie schedule time (the estimated time of arrival - the starting time) is corrpared with the calculated 
necessary time (Step 105). When the schedule time is smaller than the necessary time Qn the case where the train is 
55 too late for the necessary time at any high speed), it is judged that tiie input data are abnormal. Then, this procesang 
is completed (Step 106). The completion of this processing is displayed on the CRT 210. 

On the other hand, when tiie schedule time is larger than the necessary time, the difference between the schedule 
time and the necessary time is calculated.to compare the magnitude of the drffer^ice and the allowat^e error 6T with 
each other (Step 107). If the difference tiierebetween is smaller than tiie allowable error (formation of end condition), 
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the target speed and the running curve data at this time are stored in the running plan data file 204 (St^ 108). TTius, 

the processing for obtaining the optimal train running method is completed (Step 109). 

If the difference therebetween is more than the allowable enror 6T, the processing for changing the target speed is 

performed (Step 110). There are various methods for changing the target speed. For exanrple, there is a method in 
5 which the gradient vector having dE/dVt as the i-th componerrt is normalized (making the length thereof 1). the length 

along the i-th component, i.a. that of the target speed as the vector is decreased by 5 km/h. There is another methiod 

in which the positive and maximum partial derivative is selected and the target speed is reduced by 1 km/h. In the 

present emtxxlimerTt the latter Is employed. 

Thereafter, the sut>sequent target speed s newly set. a series of processings beginning from Step 104 are 
10 repeated. 

As another example in the method of producing a train running plan, it is supposed that the territory for providing 
the optimal train running plan exists between the two passed stations, as shown in Rg. 15. At this time, there is no 
schedule time between these stations. However, the difference between the estimated times of passing with respect to 
the two stations is regarded as tiie schedule time, tiie speed at passing tiie station at the left end (tiie initial speed in 
IS the simulation) is added, and the consumed energy and the necessary time are obtained by the simulation. 

Thus, in the method of producing a train running plan according to the first emk^odiment, since the target speed in 
each of the limit speed territories is obtained on the k)asis of the nonlinear planning for obtaining the desired train run- 
ning plan, the number of parameters used for determining tiie target speed can be decreased, and the computation 
time becomes less in comparison with the case of the optimal running metiiod for all the territories by the dynamic plan- 
20 ning or the like. 

Moreover, ttie method of producing a train running plan of the present embodiment can be applied to the running 
pattern using the arbitrary limit speed. 

The second emt>odiment of the present invention will sut>sequentiy be described. In ttie second embodiment the 
subsequent target speed is newly obtained on the basis of tiie change ratio of the consumed energy witti respect to the 
25 change of the running time when changing the target speed. Inddentally. as the system for carrying out the metiiod of 
. producing a train running plan of the second embodiment, as shown in Fig. 1 1 is used. Therefore, the description of the 
arrangement of the system is omitted here for brevity. Then, only the features of the second embodiment will hereinafter 
be described with reference to Fig. 1 6. 

In Step 101a, the data of the type of the train, the identifiers with respect to the individual stations, the starting time. 
30 the estimated time of an^ival, and the like are inputted from the keyboard 209. 

In Step 102a. on the basis of the data inputted in Step 101a, tiie line condition data with respect to ttie stations at 
which the running train is to be stopped (the line length, the station position, the limit speed infomnation, the grade inlbr- 
mation and the like) and the rolling stock characteristic data of the train of interest (the drive force characteristics, the 
deceleration characteristics, the train weight, and the like) which are stored In advance in the files 213 and 212 of the 
35 ground system 200b. respectively, are read out. The weight of the passengers in tiie train weight is calculated on the 
basis of ttie average seat-load factor (it is estimated from tiie past data). 

In Step 103a. on tiie basis of tiie data ttius read out. the initial value of tiie target speed of each of ttie limit speed 
tenritories is s^ by subtracting the marginal speed from the limit speed Vmax(n). The marginal speed is set so as not 
to exceed the limit speed during tiie fixed speed running at the maximum speed. Thus, for example, if the meirginal 
40 speed is set to 3 km/h. when the limit speed is 100 km/h, 97 km/h is set to the target speed, 
in Step 104a, the processing is performed with respect to the following two items: 

(1) The simulation of the train running is performed using tiie target speed which is set at present to calculate the 
running time T and the consumed energy E. At the same time, the combination data of the target speeds at that 

45 time, and ttie running curve data which are expressed by the relationship between ttie running position and the 
speed are temporarily stored. 

(2) The target speed of the n-tti limit speed tenitory which is set at present is expressed by V(n). Then, the target 
speed of tiiat territory is set to V(n) - eV which is obtained by reducing the target speed by the smalt value eV. to 
perform the simulation of the train running. Thus, the running time T and the consumed energy E* are calculated. 

so Then, the following approximate expression is calculated. 

p(n)^-(E-E')/{T.r) (22) 

Such a processing is canied out witti respect to n (n=1 , .... N). Incidentally, as stated in the description of ttie equa- 
ls tion (4), p(n) represents the vector component. 

Now, when the running pattern of the limit speed is ttiat as shown in Fig. 4, ttie calculation metiiod of the running 
time T and the consumed energy E is as follows. 

The running pattern of the train is supposed in such a way ttiat in the acceleration running, the acceleration by the 
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maximum accelerating force is switched over to the fixed speed running, and in the deceleration running, the decelera- 
tion by the normal maximum braking is switched over to the fixed speed running. At this time; by detern^ning the train 
running plan it means that determining the target speed in each of the limit speed territories. 

Then, if the drive force characteristics and the deceleration characteristics corresponding to the speed v are 
5 expressed by Tq(v) and Tb(v). respectively, the equation of motion of the train is given by: 

M • (dvAJt) = Fm(v) - Rr(v) - Rg(x) (23) 

where Fm(v) represents the output of the motor and is expressed by Tq(v) in the acceleration while being expressed by 
10 -7b(v) in the deceleration. Thereafter, tiie necessary time T and the consumed energy E can be obtained in such a way 
that the:time t is made to be discrete and. the above, equation is solved k>y the numerical calculation. 

Incidentally, although it is assumed that with respect to the train running, tiie acceleration running is performed in 
accordance with Tq(v) and the deceleration running is performed in accordance with Tb{v), tiie control in the fixed 
speed running is performed in the following manner. The lower limit AVI and the upper limit AVLi of tiie error between the 
15 present speed v and tiie target speed V(n) are set. Then, the control in the fixed speed running is performed in accord- 
ance with the following relationship. 

V ^ V(n) - AVI : acceleration due to Tq(^0 

20 V(n) - AVI < V < V(n) + AVu : coasting (Fm(v) = 0) 

V(n) + AVu ^ V : deceleration due to Tb(v) 

After the completion of the above processing, in^Step 105a, with the target speed being equal to the initial setting 
25 value which is set in Step 1 03a (i.a. in the.first calculation), the schedule time the estimated time of arrival - ttie start- 
ing time) is oompared. with the running time which is calculated by the simulation. 

.. . By the case where the runningtime is larger than the schedule time; it means thatin:the case where the train is too 
late for the schedule time at any high speed. In such a case, the input data are abnormal. Then, the processing for 
obtaining the optimal target speed is completed (Step 106a). Then, the completion of this processing is displayed on 
30 the CRT 210. 

In tiie case where the running time smaller than the schedule time, tiie calculation result of "the schedule time - 
the running time" and the allowable error ST are compared witti each other (Step 107a). When the difference ther^e- 
tween is smaller than the allowable en-or $T. the formation of tiie end condition is established. In this case, the combi- 
nation data of the tcuget speeds and the running curve data botii of which are temporarily stored are stored in tiie 
35 running plan data file 204 (Step 108a). Then, tiie processing for obtaining tiie optimal train running method is completed 
(Step 109a). 

On the other hand, when the difference therebetween is more than the allowable errors 5T. the processing for 
changing ttie target speed is performed (Step 11 Oa). Then, the target speed in tite territory having the positive and max- 
imum p(n) is reduced by unit quantity, e.g.. 1 km/h. After tiie subsequent target speed is newly set. a series of process- 
ed ing beginning from Step 104a are repeated. 

With respect to the method of changing the target speed in Step 1 10a, in addition to the above-mentioned metiiod. 
there are various ones which will subsequentiy be listed. 

(1) A first method is such ttiat as described on referring to the equation (4). a as the vector is used, its unit vector 
45 fi/lel is ot>talned, the target speed is moved along the direction of that unit vector by unit quantity, and in a territory 

where the target velocity exceeds the limit speed, ttiat target speed is not changed. 

(2) A second method is such that since in the above mettiod (1 ). the target speed will be increased in any territory 
having the negative p(n). a where c(n) = 0 is set to that terrtory is used. 

(3) A third metiiod is such that since in the method described in the present entediment and the atxsve methods 
50 (1) and (2). there is the possbility that the increasing of the running time oorresponding to the change of the target 

speed per unit quantity is large and thus it may not be suitably witiiin the range of the allowatde error, when 
approaching the running time set in advance, the target speed is changed to decrease the unit quantity. 

Thus, in tiie method of producing a train running plan according to the second emtxxliment the target speed in 
55 each of the limit speed tenrrtories is obtained for obtaining the train running plan, and the oomt)ination of the optimal tar- 
get ^eeds is otrtained by utilizing the ratio of ttie decreasing of the consumed wergy to the increasing of tfie running 
time in the change of the target speed. Therefore, the numb©- of parameters used for determining the target speed can 
be decreased and the computation time becomes less in comparison with the case of the metiiod of obtaining the opti- 
mal running mettiod over all the territories, such as the dynamic planning. Moreover, the present method can be applied 
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to the running pattern using the arbitrary limit speed. Further, since it is unnecessary to search for the minimum target 
speed within the limited area, the amount of calculation becomes less in comparison with the solution by the nonlinear 
planning. 

When the limit speed of the train is changed, the running plan may be corrected so that the running time becomes 
the previously determined running time. Incidentally, as the system arrangement for carrying out the method of produc- 
ing a train running plan, that shown in Rg. 11 is used. Therefore, the description of the system arrangement is omitted 
here for brevity. Then, only the features of the correction will hereinafter be described with reference to Fig. 1 7. 

In Step 101b. the data of the starting position and the ending position of the limit speed changing territory (the ter- 
ritory 5 in^Fig. 3A), arxl the change value of the limit speed are irputted. 

In Step 102b, it is judged whether or not the change of the limit speed of the limit speed changing territory influ- 
ences upon the target speed of that territory, i.e.. whether or not that target speed is allowed in the newly set limit speed. 

In the case of no influence thereupon, it is unnecessary to change the target speed, and thus the processing is 
completed (Step 103b). In the case where the influence is more or less present, a series of processing beginning from 
the Step 104b are carried out In the following manner. 

Now, the preparation for giving the description off the series of processing from Step 104b will hereinunder be per- 
formed. 

Fig. 18, Rg. 19, Rg. 20 and Rg. 21 are enlarged views of the territories 5 through 7 shown in Fig. 3A and Rg. 3B. 
Fig. 18 and Fig. 20 relate to the speed-and Fig. 19 and Rg. 21 relate to the running time density. Incidentally, it is 
assumed that the grade is zero over all the territories. The reason for setting the grade to zero was already discussed. 

Fig. 18 and Rg. 1 9 show the running plan before the reduction of the limit speed of the territory 5, and Fig. 20 and 
Fig. 21 show the running plan after the reduction of the limit speed of the ten^itory 5. In these figures, it is assumed that 
the limit speed of the territory n is VMAx,n* inverse number of V^ax j is Xmaxj, the target speed is Vn, the inverse 
number of Vn is Xn, the starting position of the territory is Xn. and the position where the switching over to the fixed 
speed running occurs is ^. In Fig. 20 and Fig. 21 , the change value is distinguished from the original value.in Fig. 18 
and Fig. 19 by putting to the original value as shown in the form of 

In Step 104b, the target speed V5' of the territory 5 where the limit speed is reduced is set to the maximum speed 
as fast as4he train runs through the territory 5. This setting is performed for the purpose of recovering the delay as soon 
as possitile to reduce the acceleration in the territory 6. Then, as the marginal speed which is set so as not to exceed 
the limit speed in the running is expressed by 6v, and then the target speed is set in accordance with the relationship of 
V5' = V,^^ 5' - 5v . If the limit speed is 1 00 km/h for example, the marginal speed Sv is set to about 3 knrt/h. 

In Step 1 05b, it is judged whether or not a limit speed territory is present between the target speed changing terri- 
tory, i.e., in this case, the territory 5 where the target speed Is changed by the change of the limit speed, and the follow- 
ing station. In the case of absence of such a territory, there is no room for correction, and it is impossible to carry out 
the scheduled operation up to the following station (Step 106b). 

In the case of presence of the limit speed territory, the processing proceeds to Step 107b. In Step 107b, it is judged 
whether or not the scheduled operation can be performed as the result of the correction of the target speed of the ter- 
ritory 6. The detail description will hereinbelow be given with respect to the judgemerrt in Step 107b. 

The running time T^ of the territories 5 through 7 shown in Rg. 18 and Rg. 19 is expressed tyy the following equa- 
tion using the equation (14) as already described. 

Ti » X4 • (xe - xg) + X4 • (X7 . xe) + G{g \Xj)) -G{g\k^)) + Xy • (Xg - ^7) (27) 

Incidentally, the position ^7 in the territory 7 where the switching over to the fixed speed running is carried out is 
obtained from the following equation by utilizing the equation (15). 

q'(V7)-q'(V4) = 47-X7 

The calculation of the running time T2 of the tenritories 5 through 7 shown in Rg. 20 and Rg. 21 will subsequently 
be performed. Now, the target speed Vg' in the territory 6 is uncertain. Then, that target speed is determined so that the 
relationship of Tg = T ^ is estatslished. 

Then, T2 is calculated while leaving V9 being an unknown quantity. Therefore. T2 is a function of Vg*. Then, while 
is and 47' are also unknown quantities, they are uniquely determined if Vg' is determined. By using the equations (14) 
and (17), the following expression is obtained. 

T2 = G(g \45')) - G(g"'{x5))+ ' (^6 - 45')+F(f "^46')) - F(f ' '(Xg)) ^^^^ 
+ >^6' •(X7 - ^6')+Q(g \47')) .G(g \x7))+ xy^ix^ - 47") 
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Incidentally, with respect to and 47*, the following relationships are established. 

q '(V5Vq*'(V4) = e5 -xs 
p''{V6Vq'{V5) = 46'-X6 

Vg' to be obtained can be calculated as the root of <MV|) = T^^y g^-T^ from the numerical calculation by utilizing 
the Newton's method or the like. 

In the case where the solution of (^(V6T=0 is presert. It Is judged in St^ 107b that it is possible to perf^^ 
scheduled operation. At this time, the target speeds V5 and Vg are corrected to V5 and Ve . respectively, to complete 
the processing (Step 108b). The resultant target speeds are displayed on the CRT 21 0, With the running curve, thecbta 
of the territories (the territories 5 through 7) where the correctbn occurs are replaced with those of the running curve 
(p(x) in Rg. 7 and q(x) in Rg. 8) when the grade is zera 

In the case of absence of the solution of <|>(V g*) = 0 , it is judged in Step 107b that it is Impossible to perform the 
scheduled operation by only tiie correction of the territories 5 and 6. At this time, in the same manner as in V5 . Vg' is 
also set to the maximum speed as fast as ttie train runs through the territory 6 (Step 109b). Then, the processing is 
returned to the Step 1 05b to perform the calculation for obtaining V/. 

In the above-mentioned algorithum. in the case where the correction is necessary for tiie fairly toward territories 
such as the territories 7_ and 8, the amount of calculation is increased so that the calculation requires much time. How- 
ever, since the target speed of tiie nearest territory is set to the maximum speed as fast as tiie frain runs through that 
territory. ,the predetermined calculation may be completed during the running through tiiat territory. Thus, tiiere is no 
problem in practical use. 

As described above, according to the method of :^oonrecting^a^ which are calculated in 

advance t)y the simulation are used to perform the calculation of the running time, thereby to perform the correction of 
.ihe~runriing pian.-::Therefore,^even. when- the change of the limit speed occurs temporarily, the running plan can be cor- 
rected in a real-time manner. 

Claims 

1. A method of produdng a train running plan for making a train to run through a predet^-mined travel territory of a 
railway, wherein tiie travel ten^itory is divided into a plurality of territories, each having a predetennined limit speed 
{Vmax,i). while maintaining a predetermined running time (Tq), tiie method comprising the steps of: 

(a) setting a target speed for each territory; 

ft)) obtaining a consumed energy (E) and a running time (T), wherein the train is simulated to run in accordance 
with the target speeds in each territory, assuming acceleration of the train t>y a maximum accelerating force 
and deceleration of the train by a maximum decelerating force: 

(c) for each territory, without altering ttie target speeds of the other territories, reducing the target speed of a 
territory and obtaining a consumed energy (E*) and a running time (T) for the changed target speed of the ter- 
ritory, arKi determining a change ratio of tiie consumed energy obtained for tiie ten-itory; and 

(d) repeating steps (b) and (c). wherein, after each step (c). the target speeds are reset on the basis of the 
change ratios, in order to determine the target speeds of all the territories in the predetermined travel territory 
in such a way that in tiie simulation, the train runs in accordance with said predetermined running time and the 
consumed energy becomes minimum. 

2. The method according to daim 1 , wherein step (c) comprises tiie steps of 

(Ci) determining the change ratio of the consumed energy for each territory as; positive for tiie case that the 
consumed energy is decreased while tiie running time s increased, wherein the target speed initially s^ in 
step (a) is assumed to be the limit speed: 

(C2) reducing the target speed of the territory, in which the resultant change ratio of the consumed energy takes 
the maximum value, by unit quantity; and 

(C3) treating the target speed reduced by unit quantity in step (C2) as a new target speed, when the mnning time 
is within a predetern^ned error range. 

3. The method according to claim 2, wherein step (C2) includes a processing of the value of the change ratio of ttie 
consumed energy obtained for each tenritory as a vector quantity to change the target speed along the direction of 
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the vector by unit quantity. 

4. The method according to daim 3. wherein said processing as a vector quantity includes a changing of a compon ent 
of the vector quantity of the territory to "0** if the component of the vector quantity is negative. 

a The method according to claim 2, wherein step (C2) includes a process of decreasing the unit quantity for chang ing 
the target speed, when the running time of the train in the ten-itory approaches the predetermined running tim&. 

6. The method according to claim 1. wherein the running time and the consumed energy of the train of the predeter- 
mined travel territory are obtained on the basis of an initial speed when the train enters into said travel territory, pre- 
determined grade information, a running resistance of the train, acceleration and deceleration characteristics of the 
train and a weight of the train. 

7. The metiiod according to claim 6. wherein the consumed energy is made to be an evaluation function (E) 
expressed by an equation of E » E(V^. .... V,^) wherein the individual target speeds are independent variables 
(Vi Vn). 

8. The method according to claim 7. wherein the running time is made to be evaluation function (T) expressed by 
equation T = T(Vi V^). 

9. The method according to claim 7 or 8, wherein a differential coefficient is obtained from an approximate expression 
of dE/dV(n) = AE/AV(n), wherein target speed V(n) Is changed in tiie range of 0 ^ V(n) ^ Vmax(n) and V{n) 
becomes an optimal target speed when E takes a mininujm value. 

PatentansprOche 

1. Ver^i^en zum Erstellen eines Zugfahrplans, um einen Zug durch ein vorbestimmtes Reisegebiet eines Eisenbahn- 
netzes fahren zu lassen. wobei das Reiseget>iet in eine Vielzahl von Get»eten unterteilt ist. die jedes eine vorbe- 
stimmte HOchstgeschwindigkeit (Vj^ax, i) besitzt wShrend eine vorbestimmte Fahrzeit (Tq) aufrechterhalten wird, 
mitden Schritten 

(a) Festiegen einer Zielgeschwindigkett fOr Jedes Gebiet* 

(b) Bestimmen der verbrauchten Energie (E) und der Fahrzeit (T). wobei simuliert wird. daB der Zug in Uber- 
einstimmung mit den Zielgeschwindigkeiten in jedem Qebiet fdhrt unter der Annahme einer Beschleunigung 
des Zuges mit einer maximalen BescNeunigungskraft und einer Abbrenrvsung des Zuges mit einer maximalen 
Abbremskraft; 

(c) fOr jedes Gebiet. Verringern der Zietgeschwindigkeit eines Gebietes. ohne die Zielgeschwindigkeiten der 
anderen Gebiete zu verandern, und Ermittein der verbrauchten Energie (E*) und der Fahrzeit (T) fQr die gean- 
derte Zielgeschwindigkeit des Gebietes und Bestimmen eines Verdnderungsquotienten der fOr das Gebiet 
ermittelten. verbrauchten Energie; und 

(d) Wiederholen der Schritte (b) und (c). wobei nach jedem Schritt (c) die Zielgeschwindigkeiten auf der GrurxJ- 
lage der Andeningsquotienten neu festgelegt werdea um die Zielgeschwindigkerten alier Gebiete in dem vor- 
beslimmten Reisegebiet auf eine solche Weise zu bestimmen. daB bei der Simulation der Zug in 
Obereinstimmung mit der vortoestimmten Fahrzeit fShrt und die verbrauchte Energie minimal wird. 

2. Verfahren gemdB Anspruch 1 . wobei der Schritt (c) die Schritte 

(Cf) Bestimmen des Anderungsquotienten der verbrauchten Energie fOr jedes Gdsiet als positiv, wenn die ver- 
brauchte Energie abnimmt. wdhrend die Fahrzeit zunimmt. wobei die anfanglich in Schritt (a) festgelegte Ziel- 
geschwindigkeit als Hdchstgeschwindigkeit angesehen wird; 

(C2) Verringern der Zielgeschwindigkeit des Gebietes. in dem der resuitierende Anderungsquotient der ver- 
brauchten Energie den maximalen Wert annimmft. um einen Einheitswert; und 

(C3) Behandein der um den Einheitswert in Schritt (cg) verringerten Zielgeschwindigkeit als neue Zielgeschwin- 
digkeit, wenn die Fahrzeit innerhalb eines vorbestimmten Fehlerbereichs liegt 

3. Verfahren gemaB Anspruch 2. wobei der Schritt (C2) ein Verarbeiten des Wertes des Anderungsquotienten der for 
ledes Gebiet ermittelten veriDrauchten Energie als VektorgrOBe. um die Zielgeschwindigkeit entlang der Richtung 
des Vekfiors um den Einheitswert zu dndern. 
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4. Verfahren gemdB Anspruch 3. wobei das Verarfoeiten einer VektorgroBe ein Andem einer Komponente der Vektor- 
grOBe des Gebietes auf 'X)" errthait. wenn die KDmponente der Vektorgr6Be negativ ist. 

5. Verfahren gemdB Anspruch 2, wobei der Schritt (Ca) einen Vorgang des Verringems des Einheitswerts zum Andern 
der Zielgeschwindigteit enthdlt, wenn die Fahrzeit des Zuges in dem Gebiet sich der vorbesUmmten Fahrzeit 
nahert. 

6. Verfahren gemdB Anspruch 1 . wobei die Fahrzeit und die verbrauchte Energle des Zuges in dem vorbestimmten 
Resegebiet auf der Gbrundlage einer Eingangsgeschwindigkeit wenn der Zug In cteis Reisegebiet einfdhrt einer 
vorbestimniten Steigungsinformation/eines:Fahnmd Abbrems- 
^characteristik des Zuges.urtd eines Gewichts des Zuges erhalten wird. 

7. Verfahren gemdB Anspruch 6. wobei die verbrauchte Energie als Auswertefunktion (E), ausgedruckt durch eine 

Gleichung E = E(V ^ V^^) . wobei die enzelnen Zielgeschwindigkeiten unabhdngige Variablen (V^ Vn) sind. 

behandelt wird. 

8. Verfahren gemdB Anspruch 7. wobei die Fahrzeit als Auswertefunktion (T), ausgedrOckt durch die Gleichung 
T = TCV , Vn) , behandelt wild. 

9. Ver^rm gemdB Anspruch 7 Oder 8. wobei der DifferentPalkoeffizient von einem Ndherungsausdruck von 
.6E/5V(n) =^E/AV(n) erhalten wlrd,.wobei die Zielgeschwindlgkeit V(n) im Bereich von 0 ^ V(n)^ V,^n) gedn- 
dert wird und V(n) die optimale Zielgeschwindigkeit ist, wenn E einen minimalen Wert annimmt 

Revendlcations 

1 . Proc6d6 de production d'un plan de circulation de train pour fatre circuler un train d travers un trongon de d^iace- 
• I -.-.ment pr6d6termin6 d^une.voiedechemin deier dans lequel le trongon de d^lacement est divisd en plisieurs tron- 

gons, ayant chacun una vitesse limite pr6d6temnin6e QJ^fio^d' tout en maintenant un temps de circulation 
pr6d6termin6 (Tq), le proc^6 comportant les stapes consistent d : 

(a) dtablir une vitesse cible pour chaque trongon ; 

(b) obtenir une dnergie consommde (E) et un temps de circulation (T), dans lesquels on simule que le train cir- 
cule confbrmdment aux vitesses cibles de chaque trongon, en supposant que le train est acc^l^d par une 
force d'accel§ration maximale et que ie train est ddcdidre par une force de ddcdl^ration maximale ; 

(c) pour chaque trongon. sans fairs varier les vitesses cibles des autres trongons. rdduire la vitesse cible du 
trongon et. obtenir une dnergie consommde (E*) et un tenps de circulation (T) pour la vitesse dbte changde 
du trongon, et ddterminer un rapport de changement de I'dnergie consommde obtenu pour le trongon ; et 

(d) rdpdter les dtapes (b) et (c). les vitesses cibles. aprds chaque dtape (c). dtant rdinrtialisdes sur la base des 
rapports de changement afin de ddterminer les vitesses dbles de tous les trongons du trongon de ddplace- 
ment prdddtermind de telle sorte que, dans la simulation, le train drcule confonmdment audit temps de circu- 
lation prdddtermind et I'dnergie consommde devient minimala 

2. Procddd selon la revendication 1 , dans lequel Tdtape (c) conporte les dtapes oonsistant d : 

(Ci) ddterminer le rapport de diangement de I'dnergie consommde pour chaque trongon comma dtant positif 
dans le cas ou I'dnergie consommde diminue alors que le terrps de circulation est augmentd, la vitesse able 
dtablle initialement d I'dtape (a) dtant supposde dtre la vitesse limite ; 

(C2) rdduire la vitesse cible du trongon. dans lequel le rapport de changement rdsultant de I'dnergie consom- 
mde prend la valeur maximale. d'une quantitd unitaire ; et 

(C3) traiter la vitesse cible rdduite d'une quantitd unitaire dans I'dtape (c^ en tant que nouvelle vitesse cible. 
forsque le temps de circulation est dans une plage d'erreurs prdddterminde. 

3. Procddd selon la revendication 2, dans lequel I'dtape {c^ comporte un traitement de la valeur du rapport de chan- 
gement de I'dnergie consommde obtenu pour cfiaque trongon en tant que grandeur vectorielle pour changer la 
vitesse cible le long de la direction du vecteur d'une quantitd unitaire. 

4. Procddd selon la revendication 3. dans lequel ledit traitement en tant que grandeur vectorielle conporte le chan- 
gement d une conposante de la grandeur vectorielle du trongon d "0* si la conposante de la grandeur vectorielle 
est ndgative. 
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5. Proc§dd selon la revendication 2, dans lequel r^tape (C2) comporte un traitement de diminution de la quantity uni- 
taire pour changer la vitesse cible, lorsque le temps de circulation du train dans le trongon approche le temps de 
circulation pr6d6termin6. 

5 6. Proc6d6 selon la revendication 1, dans lequel le temps de circulation et r^nergie consomm^e du train du tron^n 
de d^lacement prM6termin6 sont obtenus sur la base de la vitesse initiale lorsque le train p6n^re dans ledit tron- 
gon de d^placement, d*une information de d6clivit6 pr6d6termin§e, d'une r^^ance ei la circulation du train, des 
caract^ristiques d'acc^tdration et de d6c6t6ration du train et du poids du train. 

10 7. Proc6d6 selon la revendication 6, dans lequel T^nergie consomm66 est amende ^ §tre une fbnction d*6valuation 

(E) exprim^e par equation E = E(V ^ V,^) dans laquelle les vitesses cibles Individuelles sont des variables 

ind^pendantes (V^, .... Vn). 

8. ProcM^ selon la revendication 7, dans lequel le temps de circulation est amen^ k dtre une fonction d'^valuation 
15 (T) 6xprim6e par I'^uation T = T(V ^ V^). 

9. Proc6d6 selon la revendication 7 ou 8, dans lequel un coefficient diff^rentiel est obtenu & partir d'une expression 
approch6e de 6E/6V(n) = AE/AV(n) , la vitesse dble V(n) 6tant chang6e dans la plage de 0 s V(n) ^ Vmax(n) et 
V(n) 6tant une vitesse cible optimale lorsque E a une valeur minimale. 

20 



25 



30 



35 



40 



45 



SO 



55 



15 



EPO 467377 B1 

^ ^ ^- * PRIOR ART 




TRAVEL DISTANCE 



Fl G. 2 




TRAVEL DISTANCE 



16 



EPO 467377 B1 



FIG. 3A 



Q 
LU 
IjlI 
Q. 



TERRTPDRY 2. 



TERRITORY 



TERRITORY 4 TERRITORY 6 

TER RITORY 3 / TERRTTORY 5 i 




REDUCTION OF 
LIMIT SPEED 



PRESENT 
POSTION 



TERRmORY 7 



TERRrTORY 8 



TERRITORY 9 
TERRITORY 10 



POSITION X 



FIG. 3B 



>- 

tz 
UJ 

1 

o 




PRESENT 
POSITION 



17 



EP0467377B1 




18 



EPO 467377 B1 



Fl G. 5 




i 



LIMIT SPEED 



TARGET 
SPEED 



7 



RUNNING CURVE 




Fl G. 6 



POSITION X 



CO 

UJ 

Q 

111 




RUNNING TIME 
DENSITY (= l/v) 



MNVERSE NUMBER 
OF TARGET SPEED 



POSITION X 



INVERSE NUMBER 
OF LIMrr SPEED 



19 



EPO 467377 B1 



Fl G. 7 




SPEED v = P(x) 



RUNNING TIME T = F( x) = 5*f(x)dx 



RUNNING TIME DENSITY X =f(x)=l/v(x) 



POSITION X 



Fl G. 8 




SPEED v»g(x) 



RUNNING TIME DENSITY X=g(x)=l/V^(x) 



RUNNING TIME T = G( x ) =Lg( x)dx 



POSITION x 



20 



EP0467377B1 



F I G. 9 





V2 










1 




1 

i i 



XI Cl X2 ^2 X3 



FIG. 10 



XI 








X3 


















X2 


A 










B 




D 





XI C I X2 C2 



21 



EP0467 377B1 



FIG. It 



204 



RUNNING 
PROGRAM 
OAIA FILE 




203 200a 



RUNNING 


SUPPORT 


UNIT 





205 
J. 



ROLLING STOCK 
POSITION 
DETECTION UNIT 



207 



IC CARD 
READER 



206 



IC CARD 




22 



EP0 467377B1 



READING OUT THE DATA INPUTTED 
FROM THE KEYBOARD 






READING OUT THE LINE CONDITION 
DATA AND THE ROLLING STOCK 
CHARACTERISTIC DATA 







INITIALLY SETTING THE TARGET 
SPEED TO "THE LIMIT SPEED- 
THE MARGINAL SPEED " IN EACH 
OF THE LIMIT SPEED I bH W IUWES 



101 



FIG. 12 



102 



103 



OBTAINING THE NECESSARY TIME , 
THE CONSUMED ENERGY AND THE 
l=ARTIAL DERIVATIVE OF THE CON - 
SUMED ENERGY WITH RESPECT TO 
EACH OF THE TARGET SPEEDS 
IN THE SET TARGET SPEED 



104 



'WITH THE INITIAL SETTING 
VALUE, 

THE SCHEDULE TIME < 
THE NECESSARY TIME 




106 



END 

(THE INPUT DATA 
ABNORMAL > 



ARE 



NO 



'I THE SCHEDULE TIME - 

THE NECESSARY TIME I < 

THE ALLOWABLE ERROR 




NO 



REDUCING THE TARGET SPEED 
IN THE TERRITORY HAVING THE 
MAXIMUM PARTIAL DERIVATIVE BY 
Ikm/h. DO NOT CHANGE THOSE 
IN OTHER TERRITORIES 



108 



STORING THE TARGET 
SPEED AND THE RUN- 
NING CURVE DATA IN THE 
OPTIMAL TRAIN RUNNING 
METHOD DATA FILE 





109 


END 



23 




24 



EP0467 377B1 



FIG. 14 



Q 
LU 

LJ 
CL 

cn 



Vmax 



UJ 

Q 



I 

_1 

I 



rn,m 




Xn DISTANCE JC 
rptMn-l 



Cn,Mn-l 



FIG. 15 



o 

LiJ 
LU 
Q. 
CO 



Vmax 2 













1 





XO XI 



X2 



I 1 



X3 X4 X5 DISTANCE X 



25 



EP0467 377B1 



INPUTTING THE KIND OF THE TRAIN, 
THE IDENTIFIERS WITH RESPECT 
TD THE INDIVIDUAL STATIONS OF 
SUBJECT OF THE PRODUCTION, THE 



101 a 



FIG. 16 



OF INTEREST AND THE ROLLING STOCKT 

CHARACTERISTIC DATA OF THE TRAIN 
OF INTEREST 



102a 



SETTING THE INITIAL VAULE OF THE 
TARGET SPEED OF EACH OF THE 
UMIT SPEED TERRITORIES TO 
"THE UMIT SPEED -THE MARGINAL 

seeed! 



103a 



PERFORMING THE SIMULATION OF THE 
TRAIN RUNNING USING THE TARGET 
SPEED WHICH IS SET AT PRESENT 
STORING TEMPORARIiy THE TORGET 
SPEED DATA AND THE RUNNING 
CURVE DATA. PERFORMING THE SM- 
ULATION WCN THE TARGET SPEED 
I^^CH ^ THE LIMIT SPEED TER- 

VSfUJUE.AND CAlJCULAn?IG^?Jfe^!§li' 



RUNNING TIME 



104a 



< 



WTTH THE INITIALLf SET TARGET 
SPEED , 

THE RUNNISIG TIVE < THE SOEDU 



YES 



■ — ^I05a 
etne/ ho 



i06a 



END 

(THE SWT D«IA ARE 
ABNORMAL) 



107a 



iC^a 



THE SCHEDULE TIME - 

THE RUNNING TIME 

THE ALLOWABLE TIME 



NO 



REDUaNG ONLY THE TARGET 
SPEED IN THE TERRITORY H/SV- 
ING THE MAXIMUM RATD CALCU- 
LATED IN THE STEP I04a BY 1 km/h 





STORING THE TARGET 
SPEED DAIA AND THE 
RUNNING CURVE DffIA 
WHICH ARE TEMPORARILf 
STORED , IN THE OPTIMAL 
TRAIN RUNNING METHOD 
DATA FILE 


^YES 

llOa 
—/ 






I09a 

s 




END 



26 



EP0467377B1 



FIG. 17 



lOlb 
1 



INPUTTING THE DATA OF THE 
CHANGE PaNT OF THE LIMIT 
SPEED 



'IS IT NECESSARY TO CHANGE 
/THE T2XRGET SPEED OF THE 
'LIMIT SPEED CHANGING 
TERRITORY ? 



YES 




103b 



IT IS UNNECESSARti 
TO CHANGE THE 
TARGET SPEED 



CHANGING THE TARGET SPEED 
OF THE UMIT SPEED CHANGING 
TERRITORY 



IS THE LIMIT SPEED TERRI- 
TORY PRESENT BETWEEN THE 
TARGET SPEED CHANGING 

AND THE FOLLOWING 



i05b 



NO 





YES 


/can THE TARG 


ET SPEED OFN 



I06b 



IT IS IMPOSSIBLE 
TO PERFORM THE 
SCHEDULED OPER- 
ATION UP TO THE 
FOLLOWING SWnON 



I07b 



THE LIMIT SPEED TERRITTDRYV ^„ 
NEXT TO THE TARGET SPEED ^^ES 
CHANGING TERRITORY BE SET SO, 
AS TO ALLOW THE SCHEDULED 
.OPERATION ? 



108b 



SETTING OF THE 
TARGET SPEED 
IS COMPLETED 



NO 



SETTING THE TARGET SPEED OF 
THE ABOVE LIMIT SPEED TERRI- 
TORY TO THE MAXIMUM SPEED 
AS FAST AS THE TRAIN RUNS 
THERETHROUGH 



109b 

1 



27 



EP0467377B1 



F I G. 18 



VmAX.4 



LIMIT SPEED 



V4 






L X RUNNING 
\ VmAX.7 


TERRITORy 4 


■ 

TERRITORY 5 


TERRITORT 6 




;territo- 

RY 7 . 1 



X5 



X6 



X7 €7 X8 



FIG. 19 



X4 



XMAX,4 



X5 



X7 



i 



RUNNING TIME DENSITY 



"7 — I 

/ XmaX,7 



INVERSE NUMBER ! 
OF UMIT SPEED 



X6 



X7 C7 X8 



28 



EP 0467 377 B1 



FIG. 20 



Vmax,4 



V4 



UMIT SPEED 



RUNNIIM6 CURVE 



\ Vmax. 5^ 




TERRrroRY4 




TERRITORY 5 




TERRITORY 6 


TERF 


ITORY 7 j 

















ii 

X5 



3C6 ^fe 



3C7 6- X8 



FIG. 21 





X'5 


RUNNING TIME 


X7 






('DENSITY > 














XMAX,7 






'/\ Xmax, 5/ 


_\ x'6 


/ 




XMAX4 




















t 


INVERSE NUMBER 
OF LIMIT 
SPEED 


/ j Xmax.s' 

1 i 




■ 

1 





XS^'S X6 Cfe X7 ^7 



29 



